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There are many theories why people get old: programmed cell death, oxygen 
radicals, mechanical wear and tear, etc. Rose 1 suggested that methods similar to 
mathematical physics should be applied to this problem, and several Monte Carlo 
studies of biological aging have been published recently 2 . They, as well as the 
present paper, try to explain aging through population dynamics and mutations: 
Bad mutations affecting young age before reproduction reduce population growth 
much stronger than if the same mutations affect only old age after reproduction. 
This Darwinistic selection pressure thus keeps youth more free from deleterious 
mutations than old age, and thus leads to a death rate increasing with age. The 
aim of the present work is to explain the algorithm for a numerically particularly 
efficient model 3 and to see if it gives biologically reasonable results (survival rate 
decaying rapidly at old age) for parameters which require large-scale parallelization. 
This model has the advantage, compared to earlier ones 2 , that a large number of 
time intervals can easily be incorporated into the life history of an individual, in 
contrast to only two intervals (youth and adulthood). 

We start with a large number of individuals having in their genes for each "year" 
t of their lifes at most one mutation which becomes relevant in this year and the 
following years only: t = 1,2,..., 32. All mutations are deleterious (bad), and 
an individual dies if the number of mutations relevant at time t reaches a certain 
limit T. These mutations are stored as bits in a 32-bit computer word defining the 
genome of the individual. They are hereditary in the sense that at the birth of a 
child they are given on to this child, except for M randomly selected bits which are 
toggled at birth. These new mutations of rate M affect only the child in its later life, 
not the parent; and these new mutations may not only increase but also decrease 
the number of old bad mutations stored in the computer word: If the randomly 
selected bit was already set to unity before (indicating an inherited disease) it is 
now reversed to zero such that no bad mutation is stored at that age. (Sexual 
reproduction is ignored at present and could be taken into account by a genetic 
algorithm mixing two bit strings.) In each year, with probability m, one child is 
produced by each surviving individual. Survival requires that the total number of 
set bits in positions 1,2, . . .,t of the computer word is less than T. In addition, 
we let individuals with a lower number of mutations survive only with the Verhulst 
probability (1 — N/N max ) where N is the total number of invividua and N max , 
typically ten times the initial N , is a maximum number of the population due to 
food and space restrictions. 

Initially, the bits representing the genome are set randomly, half up and half 
down. The peculiar dangers of childhood are ignored here; thus babies of age t = 
should be interpreted as offspring which has overcome these initial problems. Only 
individuals older than some minimum age get children. 

This model is similar to that of Hotzel 4 in that more than two age intervals are 
taken into account. Hotzel simulated, for four aging steps, biased random walks on 
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a large four-dimensional lattice. In this sense, the present new mutations happening 
at birth correspond to random walks on a 32-dimensional hypercube with 2 32 sites. 
Our bit storage allows for a simple transfer of genetic information from parent 
to child: 32 age intervals in one operation only 5 . More biological background on 
mutations, though not for the purpose of aging, is given by Lynch and Gabriel 6 in 
connection with their Monte Carlo simulations. 

To parallelize this algorithm on a multiple-data multiple-instruction computer 
with distributed memory is quite easy, since only through the food restriction factor 
(1 — N/N max ) the different nodes of a parallel computer are coupled. The total 
number N can be easily evaluated by the global summation function gisum on the 
Intel Hypercube or Paragon computer. If the population in one of the nodes has 
died out, the loop over the predetermined number of iterations (years) still has to 
be continued since otherwise due to the synchronizing gisum all the other nodes 
would be blocked. A parallel computer with large total memory is required if, as is 
often the case, the population first decreases by several order of magnitude before 
Darwinistic selection of the fittest lets it grow again. 

In the complete Fortran program listed as Fig.l, minage is the minimum age 
for reproduction, and our variables T, N max , N , M are denoted as limit, nindm, 
npopul , mutr ; ind is the maximum number of individuals allowed because of mem- 
ory limitations. The data statement gives initial parameters. Now we explain the 
various do-loops: Loop 1 produces isolated bits set at position t = 1,2,..., 32 and 
used later for the new mutations. Loops 2 and 3 generate the random initial genes 
mgen, using multiplication with 16807 for random number generation. The main 
loop 4 counts the iteration. Loop 5 over all individuals is realized by goto 5 to 
allow for a fluctuating total number of individua. The bits (bad mutations) up to 
the current age are counted in loop 6. New mutations are introduced in loop 7 at 
birth. For parallelization we have to keep in memory the population for one node 
(loop 9) before all nodes are summed up by call gisum. At the end we print out 
the population as a function of age and the ratio between consecutive ages; the 
latter ones are the survival rates in a stationary state. 

parameter(ind=2000000) 
byte age (ind) 
real*8 dclock 

dimension mgen(ind) ,mbit (32) ,number(0 : 33) 
1 , iwork(0 : 33) ,numt (0 : 33) 

data mutr/1/, limit/ 3/, iseed/1/, max/ 1000/ 
1 , birth/0.11/, npopul/200000 / 

1 , minage/ 8/, nindm/20000000/ 

time=dclock() 

node=mynode ( i ) 

if (node . eq. 0) 

1 print 100, ind, npopul, mutr, limit, iseed, max, minage, nindm 
1 ,numnodes() , birth 
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ibm=iseed*2-l+2*node 

ibirth=2147483647 . 0*(2 . 0*birth-l . 0) 

mbit(l)=l 

number ( ) =npopul 

do 1 j=l,32 

number (j )=0 

1 if(j.gt.l) mbit(j)=ishft(mbit(j-l) ,1) 
do 2 i=l,npopul 

ici=0 

do 3 j=l,32 
ici=ishf t (ici , 1) 
ibm=16807*ibm 
3 if(ibm.lt.O) ici=ici+l 
mgen(i)=ici 

2 age(i)=0 

do 4 itime=l,max 

if (npopul . le . . or .npopul .gt . ind) goto 4 
iverhu=2147483647 . 0*(npopul*2 . 0/nindm-l .0) 
i=0 

5 i=i+l 

10 number (age ( i) ) = number(age(i) )-l 
if (age(i) .le.32) age(i)=age(i)+l 
number (age ( i) ) = number (age ( i) )+l 
k=mgen(i) 
negmut=0 
do 6 j=l,age(i) 
negmut=negmut+iand(l ,k) 

6 k=ishft(k,-l) 
ibm=ibm*16807 

if (ibm. It . iverhu. or .negmut . ge . limit . or . age(i) . eq. 33) then 

number ( age ( i ) ) =number (age(i))-l 

age (i)=age (npopul) 

mgen(i)=mgen( npopul ) 

npopul=npopul-l 

goto 10 
end if 

ibm=ibm*65539 

if (ibm. It . ibirth. and. age(i) . gt .minage) then 
npopul=npopul+l 
if (npopul . gt . ind) goto 4 
mgen(npopul)=mgen(i) 
if (mutr . eq. 0) goto 8 
do 7 k=l,mutr 
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ibm=ibm* 16807 

7 mgen(npopul)=ieor (mgen(i) ,mbit (1+ishf t (ibm,-27) ) ) 

8 age(npopul)=0 
number ( ) =number ( 0) + 1 

end if 

if (i . It .npopul) goto 5 

if (itime . le .minage+32 . or . itime . eq. (itime/minage)*minage) then 
npopt=npopul 
do 9 j=0,33 

9 numt(j)=number(j) 

call gisum(npopt , 1 , idummy) 
call gisum(numt,34,iwork) 
if (node.eq.O) 
1 print 100, itime, npopt, (numt(j) ,j=0, 32,8) ,number(l) 
endif 
4 continue 

100 format(lx,8i9,/,i9,f9.6) 

if (node.eq.O) print 102, numt 
102 format (lx, 10i7,/, lx,10i7,/, lx,14i5) 

if (node.eq.O) print 101, (numt(i)/(0.001+numt(i-l)) ,i=l,32) 

101 format(lx,13f6.3,/,lx,13f6.3,/,8f7.4) 
stop 

end 

Fig. 1. Fortran program, with parallelization commands for an Intel Paragon. 

Fig. 2 shows results for T = 3, N max = 5.12 • 10 9 , M = 1, minimum age = 
8, initial N = 0.512 • 10 9 , probability of birth 0.11, after 1000 years (not much 
changed after 10,000 years in smaller runs). We see that the population decreases 
with increasing age, first slowly, then rapidly. Thus the survival rate is first roughly 
constant near unity, and then decays rapidly towards about 1/2. 

In contrast to many other models 2,6 - but in agreement with Partridge and 
Barton 7 - the present model allows for a stable population even if all new muta- 
tions are negative (as is often assumed in the biological literature 7 ). A new mutation 
affects only the child, not the parent and the other offspring. These other children 
may thus dominate in the population and keep the average total number of mu- 
tations per individuals (number of bits set in 32-bit word) constant, avoiding the 
usual mutational meltdown 6 . Fig. 3 shows results for the case of only bad new mu- 
tations (ieor replaced by ior in the mutations) if we start with no mutations at 
all and set the birth probability as one; all other parameters are chosen as in Fig. 2 
except that the system was smaller to fit onto one workstation. Again, population 
and survival rate decay first slowly, then rapidly with increasing age. If we omit 
the mutations, on the other hand, all ages have the same survival probability (in 
general, mutations reduce the population). 
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Fig. 2. Decays of population and survival rate with age t. In a stationary state like here, the 
survival rate is the ratio of the population at age t to that at age t — 1. This simulation with 
512 million initial babies took about 12 minutes on an Intel Paragon parallel computer using 128 
nodes. (Finite size effects have been observed in the simulation with 200 thousands initial babies) . 
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Fig. 3. As fig. 2, but without mutations at the beginning and with only negative mutations later. 
This simulation took about 10 hours on an IBM Powerstation. 
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Fig. 4 shows the changes with the number of iterations ("years") for the simula- 
tions of figs. 2, using smaller systems. Variation of the above parameters of Figs. 2 
did not change drastically the qualitative behavior, provided the birth rate was 
adjusted such that the population did not decay to zero nor increased too rapidly 
towards the limit imposed by N max . 
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Fig. 4. Development of population with the number of computer interations (corresponding to 
biological generations or years) showing how a stationary state is established. 

Thus this model, similar to other bit-string models in biology 8 , gives a simple 
and efficient way to incorporate many aging intervals into simulations of a life- 
history model. 
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